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^Bounding  theorems  for  the  changes  in  elastic  energy  due  to  a  crack  in  a 
finite  body  and  crack  tip  Stress  Intensity  Factors  are  established  by  use  of  the 
variational  principles  of  the  theory  of  elasticity.  It  is  shown  that  the  change 
in  elastic  energy  due  to  a  crack  in  a  finite  body  is  larger  than  the  change  in 
energy  due  to  the  same  crack  in  an  infinite  body  with  the  same  boundary  condit¬ 
ions  when  stresses  are  prescribed  on  the  external  boundary,  and  smaller  when 
displacements  are  prescribed.  The  energy  changes  can  be  expressed  as  functions 
of  the  crack  tip  Stress  Intensity  Factors  and  for  special  loadings.  (over) 
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Abstract,  continued 


^bounds  for  single  mode  Stress  Intensity  Factors  are  obtained*  The  obtained 
inequalities  are  in  agreement  with  known  numerical  solutions  of  finite  cracked 
bodies. 
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1.  Introduction 

The  present  work  is  concerned  with  new  bounding  theorems  in  fracture 
mechanics.  The  theorems  are  established  by  the  use  of  the  variational 
Principles  of  the  Theory  of  Elasticity.  The  theorems  show  the  influence  of 
the  finite  size  of  the  cracked  body  on  the  change  in  elastic  energy  due  to 
the  crack  as  a  function  of  the  prescribed  boundary  conditions.  For  special 
loadings /  when  single  modes  occur/  it  is  shown  that  the  stress  intensity 
factors  (abbreviated  SIF)  of  a  finite  cracked  body  with  prescribed  stresses 
on  the  external  boundary  will  always  be  larger  than  the  crack  tip  SIF  of  an 
infinite  cracked  body  with  the  same  boundary  conditions  and  that  the  SIF  of 
a  finite  cracked  body  with  prescribed  displacements  on  the  external  boundary 
will  be  always  smaller  than  the  SIF  of  am  infinite  cracked  body  with  the 
same  boundary  conditions. 

2.  The  Bounding  Procedure 


The  problem  of  bounding  the  change  in  elastic  energy  due  to  a  crack  in 
a  finite  body  will  here  be  considered  in  terms  of  the  classical  extremum 
principles  of  minimum  potential  and  complementary  energies.  These  principles 
have  been  previously  used  for  bounding  the  effective  elastic  moduli  of 
heterogeneous  materials  [l]  and  recently  for  bounding  effective  moduli  of 
cracked  material  [2]  . 


on  \ 
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Finite  BQtiY  with 


S trasses. 


The  change  in  energy  due  to  a  crack  in  a  finite  body  with  prescribed 
stresses  on  the  external  boundary  will  be  bounded  by  use  of  the  principle 
of  minimum  potential  energy.  Consider  the  finite  body  of  volume  V,  con¬ 
taining  a  stress  free  crack  (fig.  la).  Traction  boundary  conditions  are 
prescribed  on  the  external  boundary  3, . 


W  ■  TI 


The  potential  energy  of  the  cracked  body  is 


a, , ,  e, . ,  dV  -  T?  u, ,  dS 


-PI  2  wijl  &ijl 

vi 


li  il 


where  and  u^  are  the  unknown  stresses,  strains  and  displacements 

in  the  body.  The  potential  energy  stored  in  the  cracked  body  is  given  by, 


UP1=°P1- 


where  Up^  is  the  potential  energy  of  the  body  of  volume  subject  to  (1) 
in  the  absence  of  the  crack,  and  is  a  function  of  0°^,  the  stresses 


and  strains  in  the  uncraoked  body, 


”?]■-$  I 


,  <n  d7 

X 
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and  5U  is  the  change  in  potential  energy  due  to  the  crack  and  according  to 
a  well  known  result,  Eshelby  [3],  can  be  expressed  as 

su  ■  i  {  'W3  (5> 

on  the 
crack 

where  [u^]  is  the  displacement  jump  accross  the  crack  faces.  From  (5)  it 
is  seen  that  only  the  displacement  in  the  immediate  vicinity  of  the  crack 
has  to  be  known  to  compute  <5U. 


The  potential  energy  functioned  is  defined  by: 


1 


dV  - 


_a  ~ 
Ti  ui 


dS 


(6) 


where  u^  is  an  admissible  displacement  field  which  must  merely  be  continu¬ 
ous  in  the  region  excluding  the  crack.  Since  in  the  present  problem  boun¬ 
dary  tractions  are  prescribed  u^  is  not  restricted  by  boundary  conditions. 
The  "stresses"  and  "strains",  ct^j  and  are  defined  by 


<4) 

(7) 

ffij  *Cijkl  *kl  (b) 


The  extremum  principle  states  that 
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»p,  “P,  ‘9> 

Thus  the  problem  of  bounding  <5u^  reduces  to  that  of  constructing  an  appro¬ 
priate  admissible  displacement  field.  To  do  this  consider  another  larger 
finite  body  of  volume  V  including  V]_  and  containing  the  same  crack,  where 


V  s  V1  +  V2  >  V1  <9) 

and  apply  the  tractions  T*  to  the  external  boundary  Sj 

Ti(S2)  *  *1  (10) 

The  displacement  solution  to  this  problem  is  u^.  The  part  of  this  dis¬ 
placement  field  contained  in  the  region  of  the  larger  body  is  an  admis¬ 
sible  displacement  field  for  the  body  of  volume  subjected  to  (1),  since 
it  is  continuous  within  the  volume  V^,  excluding  the  crack. 

\  =  ui2(jjl  V1J  HU 


The  displacement  field  u^  can  be  expressed  as  the  sum  of  the  displacement 
u°2  in  the  body  V  without  the  crack  and  the  perturbation  u!_  due  to  the 
presence  of  the  crack. 


Ui2  * 


i2 


i2 


(12) 


The  strains  and  stresses  associated  with  ui2  are 


eij2  *  eij2  +  £iJ2 


(a) 


aiJ2  s  aij2  +  aij2 


(b) 


Thus 


eij  *  eij2(in  V 


(a) 


aiJ  s  aij2(in  V 


(b) 


The  potential  energy  of  the  larger  body  is  defined  as 


P2  2  J 


V  +7 
1  2 


aij2  eij2  dV  “ 


TI  ui2  45 


and  can  be  expressed  as 


^P2  "  “  2 


J7  +v  atJ2  <SZ  '  5“2 
V1  y2 


(13) 


(14) 


(15) 


(16) 


In  order  to  bound  the  real  potential  energy  of  the  finite  body  of  volume  V. 


the  expression  0^  is  defined  by 


®P  ■  *  7  L  °w  eij2  dV  - 


(17) 
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Comparing  (16)  and  (17)  it  is  easy  to  see  that  Op  takes  the  form 


UP  *  UP2  +  2 


°ij2  eij2  dV 


(18) 


By  substituting  (11)  and  (14)  into  (6)  the  potential  energy  functional  for 
takes  the  form: 


"pi  m?L  aij2£ij2dV- 


Ti  ui2  dS 


(19) 


It  will  now  be  shown  that 


0  1  U 
P  P 


(20) 


from  which  it  then  follows  by  (8)  that 


0  ID  . 
P  Pi 


(21) 


Assuming  (20)  to  be  correct  it  can  be  written  in  the  form 


H 


7Jt+v°U2E1J2  **  JS 
V1  2  2 


-  L  T?  “12  0s  +  1  . 


°ij2  Ij  2  dV  -  T  L  lJ2tj2dV" 
2  1 

(22) 


Ti  u12  45 


or  after  rearranging 


ro  |h* 
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(aij2£ij2'H7ij2eij2)dV" 

2 


12 


dS+ 


u 


12 


dS  >  0 


(23) 


It  will  now  be  assumed  that  has  the  special  form 


T3  - 

Ti s  “j 


(24) 


both  on  S.|  and  on  S^.  This  is  known  as  homogeneous  traction  boundary  con¬ 
ditions  and  the  elasticity  stress  solutions  in  this  case  are 


ooo 
aij1  =  aiJ2  s  aij 


(25) 


for  homogeneous  elastic  bodies  of  arbitrary  shape.  The  associate  strain 
and  displacements  are 


o  o  o  ,,  o 

eij1  *  eij2  s  eij  s  ijkl  aij 


(26) 


U^  (X)  a  Ui2  (X) 


ui(,)  *  eij  xi 


It  now  follows  from  virtual  work  that 


_  ( 

a  f 

it,  *1  *  J, 

5  TI  “12^  s  -  J 

,  aiJ  dV 

v2 


(27) 


Introducing  (27)  together  with  (13)  and  (25)  into  (22)  we  obtain 
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(  "  2  aij  eIj2  +  2  aij2  £ij  +  2  °ij2  £ij 2^  dV  -  0 


(28) 


By  the  synmetry  of  the  with  respect  to  ij,icl  interchange 


o®  e»  -  (j*  cp 
ij  ij2  ij2  ij 


(29) 


and  (28)  reduces  to 


/. 


v  2  aij2  eij2  dV  ^  0 
2 


(30) 


Evidently  (30)  is  correct  as  the  integrand  is  positive  definite  since 
it  is  an  elastic  strain  energy  density.  This  proves  the  inequality  (20), 
therefore  also  the  inequality  (21).  Expressing  Op  by  (17)  and  Up^  by  (3) 
and  (4)  and  using  (25),  (21)  becomes: 


i-  t 


'V 


'ij 


<5U. 


and  therefore 


(3D 


6U1  Z  <SU2  (32) 

In  words:  whan  iJ&a  &US.  dflBaggnaflUa  -Slraas  Haundacz  conditions  AC£ 

nreaeribed  an.  Ida  asternal  boundary  at  finite  bedlea.,  Ida  energy  ndanaa 
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due.  12  ms.  acaacaas.  sL  a  crack  win  is.  lacaer  la  ihs.  saallar  fcgdau 

A  limiting  case  of  the  problem  considered  is  when  the  larger  body  of 
volume  V  becomes  infinite.  It  follows  that  the  change  in  energy  due  to  a 
crack  in  a  finite  body  subjected  to  stress  boundary  conditions  is  always 
larger  than  the  energy  change  is  an  infinite  body  subjected  to  the  same 
boundary  conditions.  Thus  Eq.  (32)  can  be  extended  to  read 

501  1  <SU2  1  5u"  (33) 

Finite  body  with  prescribed  displacements. 

The  change  in  elastic  energy  due  to  a  crack  in  a  finite  body  with 
prescribed  displacements  on  its  external  boundary  will  be  bounded  by  use  of 
the  principle  of  minimum  complementary  energy.  The  procedure  is  very  simi¬ 
lar  to  that  followed  above,  but  instead  of  am  admissible  displacement  field 
an  admissible  stress  field  has  to  be  constructed. 


Let  us  prescribe  homogeneous  displacements  on  the  external  boundaries 
S1  and  S2>  of  the  two  finite  cracked  bodies  of  volumes  V1  and  VsV^+V^  il¬ 
lustrated  in  fig.1. 


u*(S,)  =  xj 


(a) 


(34) 


Uia(S2) 


'ij  J 


Cb) 
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The  stress,  strain  and  displacement  fields  in  the  two  bodies  are 

aij1’  eijr  un  and  °ij2  £iJ2’  ui2'respectively  the  presence  of  the 
crack  and  eiji»  u°i  3,1(1  aij2’  eij2*  ui2  in  the  absence  of  tbe 

crack.  The  complementary  energies  of  the  two  bodies  are 

UC1  s  2  (  aijl£ijl  47  "  [  Til  ui  dS  "  ‘  2 
V1  s 

°C2  s  2  L  _  aij2eij2  47  '  L  Ti2  ui  48  *  "  2  L  aij2eij2  dV  +  6U2  (36) 

vl  v2  2  l 

where 


aijl£ijl  dV 


+  <5U, 


(35) 


Ti  s 


aij  "j 


(37) 


As  in  the  previous  case  the  energy  changes  <5U  are  given  by  (5). 


The  stress  solution  to  the  larger  body,  of  volume  V,  over  the  part  V1 
of  this  volume  is  an  admissible  stress  field  for  the  problem  of  the  body 
V^,  since  it  satisfies  equilibrium  everywhere  in  and  has  no  displacement 
boundary  conditions  to  satisfy,  thus 


°1)2  (1”  V 


08) 
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The  potential  energy  functional  for  the  body  V1  is  defined  as 

°C1  “if  5ij  hi  ™  *  L  Ti  ui  ^  (39) 

1  S1 

where 


T. 

1 


(40) 


and 


is  related  to  5^  by  (7b) . 


The  principle  of  minimum  complementary  energy  states  that 


°C1  *  "■ 


Cl 


(41) 


Define  Og  as 

5C  ■  -  T  f  °°U2  <i2  *  +  SD2 
From  (36)  it  follows  that  Qgcan  be  expressed  as 


(42) 


C2 


1 

2  , 


_o  o 
v  iJ2eiJ2 
v2 


dV 


H 


V  +v  aij2eij2 

VI  v2 


dV  - 


12  i 


u?  dS  + 


(43) 


+ 


1 

2 


0  e° 
ij  2eij2 


dV 
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For  boundary  conditions  (34),  (25-26)  are  true  for  elastic  homogeneous 
bodies  of  arbitrary  shape.  Expressing  the  admissible  stresses  and  the  re¬ 
lated  strains  by  (13)  and  using  the  principle  of  virtual  work  for  the  vo¬ 
lume  Vg  it  can  be  shown,  as  in  the  previous  case  that 

t>C  i  5C  (44) 

It  then  follows  that 

Uc2.0Ci  (45) 

From  (35)  and  (43)  substituted  into  (45)  we  obtain 

$U2  1  6U1  (46) 

In  words:  fflagfl  same  homogeneous  displacement  boundary  conditions 
prescribed  the  external  boundaries  of  two  finite  bodies,  the  energy 

change  due  to  the  presence  of  &  crack  will  be  smaller  in  the  small en  body. 

Again,  a  limiting  case  is  obtained  when  V  becomes  infinite:  the 
change  in  elastic  energy  due  to  a  crack  contained  in  a  finite  elastic  body 
with  prescribed  displacements  on  its  external  boundary  is  always  smaller 
than  the  change  in  elastic  energy  due  to  the  same  crack  in  an  infinite  body 
with  the  same  prescribed  displacements  at  infinity.  Eq.  (46)  extends  to 


read: 
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5U  2.  6U2  1  5U1 


3.  Bounding  si  21 L 


In  the  preceding  general  bounds  on  the  changes  in  elastic  en¬ 
ergies  due  to  the  presence  of  a  crack  in  a  finite  body  were  obtained. 
These  bounds  can  be  transformed  into  bounds  on  SIF  for  special  loadings. 


The  change  in  energy  due  to  a  Griffith  crack  imbedded  in  an  isotropic 
body  as  given  by  (5)  can  be  expressed  as  a  function  of  the  crack  tip  SIF. 
For  a  3-dimensional  body 


60  =  “  [  [Kj(a)  ♦  Kjj(a)  +  ( l-vU^U) ]  da 


When  the  cracked  body  is  2- dimensional,  only  Modes  I  and  XI  exist  and 


J  [Kj(a)  +  K^j(a)]  da 


where  a  is  half  the  crack  length. 

When  the  cracked  body  becomes  infinite  the  stress  intensity  factors  are 


,  o  1/2 

KI  a  CT22  a 


v  J/2 

KII  *  al2  a 


*111  "  a23 


.  .1/2 

•  a 
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where  <^2>  o cr°^  are  the  tensile  and  shear  stresses  applied  at  infin 
ity,  when  stress  boundary  conditions  are  prescribed,  or 


E  e°  a1/2 
E  e22  a 

-  r  r°  a1/2 
=  C  e12  a 


kIII 


r  p°  a1/2 
0  e23  a 


(a) 

(b)  (51) 

(c) 


when  displacements  of  the  kind  Uj°«  e^x^ 
dary. 


are  applied  on  the  external  boun- 


For  the  general  case  of  mixed  mode  loading,  only  the  combination 
2  2  2 

[Kj  ♦  Kjj  +  (1—  of  the  stress  intensity  factors  can  be 

bounded.  For  special  single  mode  loadings  bounds  on  single  SIF  can  be  ob¬ 
tained.  A  few  examples  will  be  considered  in  the  following. 

Rramnl m  1 ; 

a)  Stress  boundary  conditions. 

Consider  a  finite  plane  reotangle  of  dimentions  2b  and  2c  containing  a 
creek  of  length  2a.,  (fig. 2).  with  prescribed  boundary  conditions  on  the 
external  edges  of  the  form: 
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a22(x2  55  *  c)  =  a° 

an  (x1  =  ±.  b)  =  0  (52) 

ai2  (X1  s  *b)  sCT12  (x2  =  ±c)  =  0 

Since  the  loading  is  symmetric  this  is  a  mode  1  loading,  and  the  only  SIF 
existent  is  K^.  The  change  in  elastic  energy  due  to  the  crack  is 

a 

5U  -  I21  |  B^(a)  da  (53) 

o 


(56) 
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_  v  o  1/2 
Kj  Iff  a 


boundary  conditions 


Consider  again  the  rectangle  of  Fig  2.  with  prescribed  displacement 
boundary  conditions  on  the  external  edges  of  the  rectangle  of  the  fora: 


u2  (x2  S  ±  C)  s  ±  vo 

U1  (x2  s  =  0 

u2  C*1  *  i  b)  *  u1  (Xj  s  i  b)  s  0 


This  deformation  is  symmetric  so  only  the  mode  I  SIF  will  enter  the 

energy  change  expression.  We  can  express  the  prescribed  displacement  v  as 

o 


o 

V0  =  e  c 


The  stress  intensity  factor  for  an  infinite  body  with  prescribed  displace¬ 
ments  at  infinity  according  to  (51a),  is 


t"  ,  E£°.1/2 
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When  the  changes  in  elastic  energy  for  the  finite  and  infinite  rectangle 
are  expressed  as  a  function  of  the  mode  I  SIF  and  inequality  (47)  is  used 
the  following  bound  on  K^(a)  is  obtained: 


Kj(a)  i  K*  (a) 


(60) 


or 


Kj(a)  i. 


Ev 


1/2 


(61) 


Isida  [7]  calculated,  numerically,  the  mode  I  SIF  for  a  finite  rectan¬ 
gle  with  various  boundary  conditions.  For  boundary  conditions  (51)  his  re¬ 
sults  satisfy  (57),  i.e.  the  SIF  are  always  larger  than  the  one  for  an  in¬ 
finite  cracked  body. 


A  second  problem  considered  by  Isida  was  that  of  the  finite  rectangle 
with  the  following  boundary  conditions 


u2  (x2  s  *  c)  s  ±  vo  (a) 

Uj  (x2  a  ±  c)  a  0  (b) 

a11  (x1  »tb)  »  d12  (x1  a  X  b)  a  0  (c) 


(62) 


It  should  be  pointed  out  that  the  results  for  this  case  cannot  be 


bounded  by  the  theorems  developed  in  the  present  work.  Since  the  boundary 

conditions  are  mixed  [  (62c)  are  traction  boundary  conditions].  However 

Evo  1/2 

these  results  also  satisfy  (61),  i.e.  is  always  smaller  than  a . 

It  should  be  emphasized  that  there  are  no  similar  inequalities  for  the 
SIF  of  a  finite  body  with  mixed  boundary  conditions.  It  may  be  either 
larger  or  smaller  than  the  SIF  for  the  corresponding  infinite  body.  Isida 
[4]  also  considered  the  case  of  the  rectangle  with  the  following  loading 

u2  (x2  z  ±0)  =  ±  vQ 


a12  (x,  =  ±  b)  =  a12(x2  =  ±.  c)  =0 

(63) 

(x,  =  ±  b)  *  0 

For  this  case  the  results  obtained  where  sometimes  larger  than  and 
sometimes  smaller,  depending  on  the  ratios  a/b  and  c/b. 

Examnla  2: 

A  second  example  considered  will  be  that  of  a  finite  rectangle  of  di- 
mentions  b  and  2h  containing  a  cracx  of  length  a,  (fig. 3). 

The  corresponding  infinite  body  for  this  problem  will  be  that  of  a 


half  plane  containing  a  crack.  The  mode  I  SIF  for  the  infinite  body  is 


1.12  a°  a1/2 


(64) 


when  a  stress 


o 

a 


perpendicular  to  the  crack  is  prescribed  at  infinity  and 


1.12  E  e°  a1/2 


(65) 


when  a  displacement  field 


(66) 


is  prescribed  at  infinity. 

Various  solutions  to  the  problem  of  the  finite  rectangle  with  boundary 
conditions 


0^2  *  ±.  h)  8  a° 

CJ12  (x,  s  0,b)  =  a12  (x2  =  4.  h)  =  0  (67) 

"it  (ii  *  °-b>  *  0 

have  been  summarized  in  [8]  and  the  K^(a)  obtained  are  always  larger  than 
1.12  <3°  a*/2  as  predicted  by  eq.  (56). 


finite  bar 


Also  the  mode  III  loading  with  prescribed  stresses  for  a 
with  an  edge  crack  was  solved  numerically,  [5],  and  the  results  are  also  as 

oo 

expected:  the  obtained  (a)  are  larger  than  . 

4.  fiaofiiiiaiflna 

By  using  the  variational  principles  of  the  theory  of  elasticity  it  has 
been  shown  that  the  change  in  energy  due  to  a  crack  in  a  finite  body  incre¬ 
ases  as  the  dimensions  of  the  body  decrease  and  is  always  larger  than  the 
change  in  energy  due  to  a  crack  in  an  infinite  body,  when  stress  boundary 
conditions  are  prescribed.  When  displacement  boundary  conditions  are  pres¬ 
cribed  the  inverse  is  true.  As  the  dimensions  of  the  body  decrease,  the 
change  in  elastic  energy  also  decreases  and  the  change  in  energy  due  to  a 
crack  in  an  infinite  body  is  always  larger  than  in  a  finite  one. 

The  changes  in  elastic  energy  can  be  expressed  as  a  function  of  the 
crack  tip  SIF  and  for  special  loadings  the  SIF  can  be  shown  to  be  always 
larger  (for  prescribed  stresses)  or  Smaller  (for  prescribed  displacements) 
than  the  SIF  of  a  crack  in  an  infinite  body. 

The  changes  in  elastic  energy  can  be  related  to  the  crack  criticality 
criterion.  The  decrease  in  the  size  of  a  finite  body  will  precipitate  fai¬ 
lure  when  stresses  are  applied  and  postpone  it  when  displacements  are  ap¬ 
plied.  When  alternating  stresses  are  applied  to  a  finite  fatigue  specimen 
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and  a  crack  is  growing  during  the  fatigue  process,  the  crack  tip  SIF  incre¬ 
ases  with  the  dimensions  of  the  crack  and  so  does  the  crack  criticality. 

Existing  solutions  to  crack  problems  in  finite  bodies  have  been  shown 
to  be  in  agreement  with  the  theorems  developed  in  the  present  work. 
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figure  fianUana 

Fig.  1:  Finite  cracked  bodies  of  arbitrary  size 

Fig.  2:  Finite  rectangle  with  central  crack 

Fig.  3:  Finite  rectangle  with  edge  crack 


